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Evidence for the neutron-rich hypernucleus A H is presented from the FINUDA experiment at 
DA«I>NE, Frascati, studying (ir + , 7r _ ) pairs in coincidence from the K^ op + 6 Li — > A H+ir + production 
reaction followed by A H — > 6 He + n~ weak decay. The production rate of undergoing this two- 
body 7T~ decay is determined to be (2.9±2.0) • W~ 6 / K^ op . Its binding energy, evaluated jointly from 
production and decay, is Ba(aH) = (4.0± 1.1) MeV with respect to 5 H + A. A systematic difference 
of (0.98 ± 0.74) MeV between Ba values derived separately from decay and from production is 
tentatively assigned to the A H 0^ s . — > 1 + excitation. 



PACS numbers: 21.80.+a, 25.80.Nv, 21.10.Gv 

Introduction. The existence and observability of 
neutron-rich A hypernuclei was discussed back in 1963 
by Dalitz and Levi-Setti pQ who predicted the stability 
of ^ H consisting of four neutrons, one proton and one 
A hyperon. Accordingly, the A hyperon stabilizes the 
core nucleus 5 H which is a broad resonance 1.7 MeV 
above 3 H + 2n [2 . To be stable, must lie also be- 
low + 2n which provides the lowest particle stability 
threshold. This motivates a + 2n two-neutron halo 
cluster structure for ^H, with binding energy and excita- 
tion spectrum that might deviate substantially from the 
extrapolation practised in Ref. pQ ■ Specifically, the study 
of °H and of heavier neutron-rich hypernuclei that go 
appreciably beyond the neutron stability drip line in nu- 
clear systems could place valuable constraints on the size 



of coherent AN — SA mixing in dense strange neutron- 
rich matter [3] . This mixing provides a robust mechanism 
for generating three-body ANN interactions, with imme- 
diate impact on the stiffness/softness of the equation of 
state for hyperons in neutron-star matter, as reviewed 
recently in Ref. [3]. 

In this Letter we report on a study of in the double 
charge exchange reaction at rest 

A; top + 6 Li -> \R + 7T+ (p,+ ~ 252 MeV/c) (1) 

based on analyzing the total data sample of the FIN- 
UDA experiment during 2003-2007 and corresponding 
to a total integrated luminosity of 1156 pb _1 . A first 
analysis of partial data, corresponding to an integrated 
luminosity 190 pb~ 4 , gave only an upper limit for (nj): 
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(2.5 ± 0.4 stat t Q A lsyst ) ■ 10- 5 /i^op 0. Although the 
statistics collected on 6 Li targets is improved by a factor 
five with respect to the run of the earlier search, the inclu- 
sive 7r + spectra do not show any clear peak attributable 
to near p n + ~ 252 MeV/c. Exploiting the increased 
statistics, the essential idea of the present analysis was to 
reduce the overwhelming background events in reaction 
by requiring coincidence with ir~ mesons from the 
two-body weak decay 

6 A R^ e Rc + TT- (p„- ~ 134 MeV/c), (2) 

with a branching ratio of about 50% considering the value 
measured for — > 4 He + ir~ [B] . The analysis described 
below yielded three distinct candidate events which 
give evidence for a particle-stable with some indica- 
tion of its excitation spectrum. The deduced binding 
energy does not confirm the large effects conjectured in 
Ref. 0. 

Data analysis. We first recall the experimental fea- 
tures relevant to the present analysis. For 7r + with mo- 
mentum ~ 250 MeV/c the resolution of the tracker was 
determined by means of the peak due to monochromatic 
(236.5 MeV/c) fj+ from decay and is a p = (1.1±0.1) 
MeV/c [7J, the precision on the absolute momentum cal- 
ibration is better than 0.12 MeV/c for the 6 Li targets, 
which corresponds to a systematic deviation on the ki- 
netic energy 0r sys t(7r + ) = 0.1 McV. For 7r~ with momen- 
tum ^130 MeV/c the resolution and absolute calibra- 
tion were evaluated from the peak due to monochromatic 
(132.8 MeV/c) ir~ coming from the two-body weak de- 
cay of iH, produced as hyperfragment with a formation 
probability about 10~ 3 — 10~ 2 per stopped K~ [BJ. A 
resolution a v — (1.2 ± 0.1) MeV/c and precision of 0.2 
MeV/c were found, corresponding to a systematic devi- 
ation of the kinetic energy o-Tsyst^ - ) = 0.14 MeV. 

Since the stopping time of in metallic Li is shorter 
than its lifetime, both production (JlJ and decay ^ oc- 
cur at rest, and a straightforward algebra leads to the 
following expression for T sum = T(n + ) + T(n~): 

T sum = M{K~) + M(p) - M(n) - 2M(tt) 

- B( 6 Li) + 5( 6 He) -T( 6 He) -T(^H), (3) 

in which M stands for known masses, B for known nu- 
clear binding energies, and T for kinetic energies. The 
evaluation of T(^H) using momentum and energy conser- 
vation depends explicitly on the knowledge of -Ba(^H), 
whereas T( 6 He) depends only implicitly on Ba(%H) 
through the momentum . 

We assume B\(\H) = 5 MeV, the average of 4.2 and 
5.8 MeV predicted in Refs. [TJ|3], respectively, with re- 
spect to 5 H + A. This choice is not critical, since T sum 
varies merely by 50 keV upon varying B\(jjl) by 1 MeV, 
negligibly low with respect to the experimental energy 
resolutions <7 T (7r + ) = 0.96 MeV and a T {^~) = 0.84 MeV 



for p^+ k, 250 MeV/c and p n - « 130 MeV/c. There- 
fore, the FINUDA energy resolution for a tt^ 1 pair in co- 
incidence is <tt = 1-28 MeV. Evaluating the r.h.s. of 
Eq. ^ one obtains T sum = 203 ± 1.3 MeV for can- 
didate events. In practice we have focused on events 
in the interval T sum = 203 ± 1 MeV, corresponding to 
only ~ 77% of the FINUDA total energy resolution; this 
value was chosen as a compromise between seeking to re- 
duce contamination from background reactions discussed 
in more detail below, and maintaining reasonable statis- 
tics, which resulted in a somewhat narrower interval than 
the experimental resolution. The raw spectrum of T sum 
for 7r^~ pair coincidence events is shown in Fig. [l] where 
events satisfying T sum = 203 ± 1 MeV are indicated by a 
vertical (red) bar. 




T(it*) + T(it ) (MeV) 

FIG. 1. (color online). Distribution of raw total kinetic energy 
T sum = T{tt + ) + T(n~) for pair coincidence events from 
6 Li targets. The vertical (red) bar represents the cut T sum = 
202 — 204 MeV. The dashed (blue) histogram is a quasi-free 
simulation of K~ op + 6 Li £+ + 4 He + n + 7r _ ; £+ ^ n + 7T+ 
background and the dotted (violet) histogram is a four-body 
phase space simulation of the same background, their best fit 
to the data is shown by the solid (black) histogram, see text. 

Figure [2] (left) shows a 2-d plot in the p v ± plane of 
coincidence events selected in the band T sum = 202 — 204 
MeV. The distribution falls to zero at p n + ~ 245 MeV/c 
and higher, and at p w - ~ 145 MeV/c and lower. This is 
close to where events are expected. Thus, to search 
for particle-stable events below its (^H + 2n) lowest 
threshold, using the two-body kinematics of Eqs. and 
([2]), a further requirement of p v + > 251.9 MeV/c and 
Pn- < 135.6 MeV/c is necessary. In the final analysis we 
selected p„+ = (250-255) MeV/c andp T - = (130-137) 
MeV/c, thus covering a mass range from the (A + 
3 H + 2n) threshold, about 2 MeV in the continuum, 
down to a bound somewhat stronger than predicted 
by Akaishi et al. [3]- This does not completely exclude 
eventual contributions from the production and decay of 
(^H + 2n) as discussed below. 

Results. Out of a total number of ~ 2.7 ■ 10 7 K~ 
detected at stop in the 6 Li targets, we found three events 
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FIG. 2. (color online), n momentum vs n momentum for Li target events with T sum — 202 — 204 MeV (l.h.s.) and with 
Tsum = 200 — 206 MeV (r.h.s.). The shaded (red) rectangle on each side consists of a subset of events with p n + — 250 — 255 
MeV/c and p„- = 130 - 137 MeV/c. 



that satisfy the final requirements: T sum = 202 — 204 
MeV, p„+ = 250 - 255 MeV/c and p„- = 130 - 137 
MeV/c, as shown within the shaded (red) rectangle in 
Fig. 2 left. Different choices of T sum interval widths (2 — 6 
MeV) and position (center in 202-204 MeV), and of p„± 
interval widths (5 — 10 and 8—15 MeV/c respectively) 
with fixed limits at 250 and 137 MeV/c respectively to 
exclude the unbound region, do not affect the population 
of this selected rectangle. For example, no new candidate 
events appear in the shaded rectangle upon extending the 
cut T sum = 202-204 MeV in the l.h.s. of Fig. 2 to T sum = 
200 - 206 MeV in the r.h.s. of the figure. A similar 
stability is not observed in the opposite corner of Fig. 2 
where, on top of the events already there on the l.h.s., six 
additional events appear on the r.h.s. upon extending the 
T sum cut of the l.h.s. Quantitatively, fitting the projected 
7r distributions of Fig. 2 left by gaussians, an excess 
of three events in both p n ± distributions is invariably 
found, corresponding to the shaded (red) rectangle. The 
probability for the three events to belong to the fitted 
gaussian distribution is less than 0.5% in both cases. This 
rules out systematic errors associated with the present 
analysis selection. 

The three candidate events are listed in Table [i] to- 
gether with nuclear mass values derived separately from 
production (JlJ and from decay These mass values 
yield a mean value M QH) = 5801.4 ± 1.1 MeV, jointly 
from production and decay, where the error reflects the 
spread of the average mass values for the three events, 
and is larger than the 0.96 MeV and 0.84 MeV measure- 
ment uncertainties in production and decay, respectively, 
for each of the three events. We note that the mass value 
inferred from the third event by averaging on produc- 
tion and decay is about 2a from the mean mass value, 



TABLE I. Summed kinetic energy T sum = T(tt+) + T(ir~), 
pion momenta p w ± , and mass values inferred for the three 
aH candidate events from production |l]) and decay The 
mean mass value is M(%B) = 5801.4 ± 1.1 MeV, see text. 



Tsum P„+ Pn- M(iH) prod . M(iH) d ecay 

(MeV) (MeV/c) (MeV/c) (MeV) (MeV) 

202.6±1.3 251.3±1.1 135.1±1.2 5802.33±0.96 5801.41±0.84 
202.7±1.3 250.1±1.1 136.9±1.2 5803.45±0.96 5802.73±0.84 
202.1±1.3 253.8±1.1 131.2±1.2 5799.97±0.96 5798.66±0.84 



an observation that could indicate some irregularity in 
the reconstruction of the third event. To regain confi- 
dence, each one of the three events was checked visually 
for irregularities, but none was found. 

Furthermore, we note from Table [I] that the mass val- 
ues associated with production are systematically higher 
than those associated with decay, by 0.98 ± 0.74 MeV 
recalling the 1.28 MeV uncertainty for T sum from which 
each of these mass differences is directly determined. Un- 
like the mean mass value, the spread of the produc- 
tion vs decay mass differences is well within la. These 
mass differences are likely to be connected to the excita- 
tion spectrum of as discussed below. 

Background estimate and production rate A 

complete simulation was performed of K^ op absorption 
reactions on single nucleons, as well as on correlated few- 
nucleon clusters, that lead to the formation and decay of 
A and £ hyperons. Full details will be given elsewhere, 
here it is sufficient to focus on two chains of reactions 
likely to produce coincidences overlapping with those 
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selected to satisfy A H production |TJ and decay 
(i) S + production 

K~ top + 6 Li -> £+ + 4 He + n + tt" , 



(4) 



where p w - < 190 MeV/c, followed by E + decay in flight 

£+ -> n + vr+ [p^+ < 282 MeV/c]. (5) 

The S + production was treated in the quasi-free ap- 
proach, following the analysis of the FINUDA experi- 
ment observing Jy^n^ pairs [5]. This simulation is shown 
in Fig. [T] normalized to the experimental distribution 
area. It provides too sharp a decrease in the 200-210 
MeV region. To have a satisfactory description a contri- 
bution (~ 25%) from a pure 4-body phase space mech- 
anism was added and a fair agreement was obtained 
(X 2 = 40/39) in the 180-220 MeV range. The simulated 
background spectra reproduce reasonably the projected 
distributions of momentum too, showing in particular 
only a small contribution to the signal region, evaluated 
to be 0.16 ± 0.07 expected events (BGD1). 



(ii) A H production 

Kto V + 6 Li 



1H 



2n + 7T+ 



(6) 



where p 7T + < 252 MeV/c, with A H decay at rest 



lH — >■ 4 He 



fa- ~ 132.8 MeV/c]. (7) 



The 7T~ momentum in this decay is close to p n - ~ 
134 MeV/c from the two-body decay of ^H, evaluated 
assuming Ba( a H) = 5 MeV as discussed above. A value 
of 0.04 ± 0.01 expected events for the ([6])-@ reaction 
chain, negligible when compared to BGD1, was obtained 
under most pessimistic assumptions for the various terms 
of the calculation. 

All other reaction chains that could produce ^ coin- 
cidences within the described selection ranges were ruled 
out by the selections applied. Turning to potential instru- 
mental backgrounds, we note that these could result from 
fake tracks, misidcntificd as true events by the track re- 
construction algorithms. To this end we considered, with 
the same cuts, events coming from different nuclear tar- 
gets used in the same runs ( 7 Li, 9 Be, 13 C, 16 0). We found 
one event coming from 9 Be. Furthermore, we considered 
events relative to the 6 Li targets, selected with a value of 
Tsum = 193 — 199 MeV, so as to search for neutron-rich 
hypernuclei produced on the other targets. No event was 
found. We evaluate 0.27±0.27 expected fake events from 
6 Li, due to instrumental background (BGD2). 

To recap, the estimated number of events due to phys- 
ical and instrumental backgrounds feeding through the 
selection criteria are 0.16 ±0.07 (BGD1) and 0.27 ±0.27 
(BGD2), giving a total background of 0.43 ± 0.28 ex- 
pected events. Thus, using Poisson distribution, the 
three A H-assigned events do not arise from background 
to a confidence level of 99%. The statistical significance 



of the result is S=7.1 considering only the physical back- 
ground, S=3.9 considering both physical and instrumen- 
tal backgrounds. 

Given the above background estimates, plus efficiency, 
target purity and cut estimates, it is possible to evaluate 
the product R(n + ) • BR(tt~), where R(tt + ) is the A H 
production rate per K^ op in reaction bin and BR(7r _ ) 
the branching ratio for the two-body n~ aecay pi): 



R{tt + ) ■ BR(tt-) = (2.9 ± 2.0) • 10" 6 /^ 



stop ' 



(8) 



Details will be given in a separate report. Assuming 
BR(7r") = 49%, as for the analogous A H 4 He + ir~ 
decay 0, we find R{tt+) = (5.9 ±4.0) • 10- 6 /i^ op , fully 
consistent with the previous FINUDA upper limit [5]. 

Discussion and Conclusion. Table [T] yields a mean 
value £ A ( A H) = 4.0 ± 1.1 MeV with respect to 5 H + A, as 
shown in Fig. [3j in good agreement with the estimate 4.2 
MeV [1 and close to B A ( A He) = 4.18 ± 0.10 MeV (with 
respect to 5 He + A) for the other known A ~ 6 hypernu- 
cleus [9], but considerably short of Akaishi's prediction 
Bf( 6 A R) = 5.8 MeV g]. This indicates that coherent 
AV — £ N mixing in the s-shell hypernucleus A H [TU] be- 
comes rather ineffective for the excess p shell neutrons in 
^H. Indeed, recent shell-model calculations by Millener 
indicate that AN — UN mixing contributions to B\ and 
to doublet spin splittings in the p shell are rather small, 
about (10 ± 5)% of their contribution in [TT] , 
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FIG. 3. (color online). A H mass (r.h.s.) from three |H candi- 
date events, as related to several particle stability thresholds 
and theoretical predictions (l.h.s.). 



Next, we ask whether the three events that give evi- 
dence for a particle-stable A H provide additional infor- 
mation on its excitation spectrum which is expected to 
consist of a + g.s. and 1 + excited state as in A H (1.04 
MeV), and a 2 + excited state as for the p-shell dineutron 
system in 6 He (1.80 MeV). In fact, it is a h (1 + ) that is 
likely to be produced in reaction ([T]) simply because Pauli 
spin is conserved in production at rest, and the Pauli spin 
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of 6 Li is S = 1 to better than 98% [TJ. The weak decay 
([2]), however, occurs from ^H(0 + ) g.s. since the (unseen) 
7 transition 1 + — ^ + is about three orders of magnitude 
faster than weak decay. Indeed, the production vs decay 
mass difference 0.98 ±0.74 MeV extracted from the three 
events listed in Table [i] is comparable to the under- 
lying 1.04 MeV 1 + excitation in but, again, smaller 
than the 2.4 MeV predicted by Akaishi et al. 3J. If this 
is the case, then the value for the g.s. would be larger 
by 0.5 MeV than that determined above, amounting to 
£?A(^H g s .) = 4.5 ±1.2 MeV. This scenario requires fur- 
ther exploration, experimental as well as theoretical. 

In conclusion, we have presented the first evidence 
for heavy hyper-hydrogen ^H, based on detecting three 
events shown to be clear of instrumental and/or physical 
backgrounds. The derived binding energy of limits 
the strength of the coherent AN — E7V mixing effect pre- 
dicted in neutron-rich strange matter [3] and together 
with the conjectured + — 1 + doublet splitting it places 
a limit on this mixing that could orient further explo- 
rations of other neutron-rich hypernuclei. A search of 
and ^Li in the (ir~,K + ) reaction at 1.2 GeV/c on 
6 Li and 10 B, respectively, is scheduled in the near future 
at J-PARC. 
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